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Abstract 

In this paper, we studied the existence theorems and techniques for finding the 
solutions of a system of nonlinear set valued variational inclusions in Hilbert spaces. To 
overcome the difficulties, due to the presence of a proper convex lower semicontinuous 
function (p and a mapping g which appeared in the considered problems, we have 
used the resolvent operator technique to suggest an iterative algorithm to compute 
approximate solutions of the system of nonlinear set valued variational inclusions. The 
convergence of the iterative sequences generated by algorithm is also proved. 

AMS Mathematics subject classification: 49J40; 47H06 
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Introduction 

It is well known that variational inequality theory and complementarity problems are 
very powerful tools of current mathematical technology. In recent years, the classical 
variational inequality and complementarity problems have been extended and general- 
ized to study a large variety of problems arising in economics, control problems, contact 
problems, mechanics, transportation, equilibrium problems, optimization theory, non- 
linear programming, transportation equilibrium and engineering sciences, see (Aubin 
1982; Baiocchi and Capelo 1984; Chang 1984; Giannessi and Maugeri 1995). Hassouni and 
Moudafi 2001 introduced and studied a class of mixed type variational inequalities with 
single valued mappings which was called variational inclusions. Since many authors have 
obtained important extension generalizations of the results in (Hassouni and Moudafi 
2001) from various directions, see (Agarwal et al. 2011; Fang et al. 2005; Kassay and 
Kolumban 2000; Petrot 2010). Verma 1999; 2001a introduced and studied some system 
of variational inequalities with iterative algorithms to compute approximate solutions in 
Hilbert spaces. 

Inspired and motivated by the research work going on this field, in this works, the 
methods for finding the common solutions of a system of nonlinear set valued variational 
inclusions involving different nonlinear operators and fixed point problem are considered 
and studied, via proximal method in the framework of Hilbert spaces. 

Since the problems of a system of a nonlinear set valued variational inequalities and 
fixed point are both important, the results present in this paper are useful and can be 
viewed as an improvement and extension of the previously known results appearing in 
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literature, which are improves the results of Chang et al. 2007 and also extends the results 
of Verma 2001b; 2002, Ahmad and Salahuddin 2012, Ding and Luo 2000, Inchan and 
Petrot 2011, Kim and Kim 2004, Kim and Hu 2008, Nie et al. 2003 and Suantai and Petrot 
2011, etc. 

Let H hea real Hilbert space whose inner product and norm are denoted by (■, ■> and 
II ■ II respectively and K be a nonempty closed convex subset of H. Let CB(H) be the family 
of all nonempty closed convex and bounded sets in H and (p : H (— oo, +oo) be a 
proper convex lower semicontinuous function on H. Let Nf.HxH^Hhea nonlinear 
function, gi : K ^ H he a nonlinear operator, Ai,Bi : K CB(H) be the nonlinear 
set valued mappings and let r, be a fixed positive real number for each i = 1, 2, 3. Set 
S = {Ni,N2,N:i],^ = {Ai,A2,A3],^ = {81,82,83}, A = {^1,^2,^3}- The system of 
nonlinear set valued variational inclusions involving three different nonlinear operators is 
defined as follows: 

Find {x*,y*,z*) € H x H x H, ul e A^ix*), e Bsix*), m| e A2(z*), e Btiz*), u\ e 
Ai{y*),vl e Biiy*), such that 

(riATi {u\,vl) +gi{x*) -gi{f),gi{x) -giix*)) - n<t>{gi{x*)) + n^giix)) > 0, gi(x) e K, 

{r2N2 (ul.vl) +g2(y*) -g2(z*),g2(x)-g2(y*)) " r24>(g2(y*)) + r2Ct>(g2(x)) > 0, g2(x) e K, 
{rsNs (ul.vl) +g3{z*) - g3{x*},g3{x} -gsiz*)) - rscpigsiz*)) + rscpigsix)) > 0, gsix) G K. 

(1) 

We denote the set of all solutions (x* , y* , z* , u^, vl, U2, v^, u^, vl) of problem (1) by 
SNSVVID(S,2l,«B,A,/<0. 
We first recall some basic concepts and well known results. 

Definition 1. A mapping g : H ^ H is said to be 

(i) monotone, if 

igix) - giy), x-y) >0 Wx,y €H; 

(ii) strictly monotone, ifg is monotone and 
{g{x) — giy), x — y) =0 if and only ifx = y; 

(Hi) v-strongly monotone, if there exists a constant u >0 such that 

(g{x) - g{y),x- y) > v\\x-y\\^, Vx,y € H; 
(iv) Lipschitz continuous, if there exists a constant v > Q such that 

\\g(x) -g{y)\\ < v\\x-y\\, Vx.y e H. 

Definition 2. A set valued mapping A : H ^ 2^ is said to be v-strongly monotone, if 
there exists a constant v > Q such that 

(wi — W2,x — y) > u\\x—y\\^, Wx,y €H,wi e A(x),W2 e Ay. 

Definition 3. A set valued mapping A : H — CB(H) is said to be x -Lipschitz continuous 
if there exists a constant t > 0 such that 

H(Ax,Ay) < t||;«; - jll, Wx,y e H, 

where Hi-, •) is the Hausdorff metric on CB(H). 
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Definition 4. (Brezis 1973) 

If M is maximal monotone operator on H then for any A, > 0 the resolvent operator 
associated with M is defined by 

}m{x) = (/ + XM)-'^(x), Wx € H. 

It is well know that a monotone operator is maximal iff its resolvent operator is 
defined every where. Furthermore the resolvent operator is single valued and nonexpan- 
sive. In particular the subdifferential dcp of a proper convex lower semicontinuous function 
(p -.H ^ (—00, +oo) is a maximal monotone operator. 

Lemma 1. (Brezis 1973) The points u,z &H satisfies the inequality 

{u — z,x — u) + X<j){x) — X(p(u) > 0, V^v e H, 

if and only if 

u=J^{z), 

where — {I -\- Xd(j))~^ is a resolvent operator and X > Oisa constant. 
For any x,y e H, is nonexpansive, i.e., 

II/^W-^WII < \\x-y\\, ^x,y€H. 

Assume that g : H ^ H is a surjective mapping and from Lemma 1 and (1) we have the 
following proximal point problem: 

' giix*) =Jl'[giiy*)-nN,{ul,vl)\, 

g2(y*) = [g2(Z*) - r2N2 {u\, V*)] , (2) 
^ ^3(z*) = \gzix*) - r^N^ {ul, vl)] , 
provided K c gi{H) for each i = 1,2, 3. 

Lemma 2. (Weng 1991) 
Let {a„], {b„] and {c„} be three sequences of nonnegative real numbers such that 

an+\ < (1 - t„)a„ +b„ + c„ V« > no, 

where no is a nonnegative integer, {t„} is a sequence in (0, 1) with Yl'r^o ^« — 
lim„_>.oo b„ = Oit„) and Yln^o < + oo. Than a„ ^ Oasn ^ +oo. 

Definition 5. Let A,B : H ^ 2^ be set valued mappings and N : H x H ^ H be a 

nonlinear mapping. 

(i) N is said to be A -strongly monotone with respect to the first argument, if there 
exists a constant v > 0 such that for all x,y e H 

[N{ui, w) — N{U2, w),x —y)>v\\x— y\\'^ Vmi e A{x), u% e A(y), w € H; 

(a) N is said to be B-relaxed monotone with respect to the second argument, if there 
exists a constant^ > 0 such that for all x, y € H,vi e B{x),V2 € B(y) 

{N{u,vy) - N{u,V2),x - y) > -^Wx-yW"^, Vu e H. 

Main results 

We begin with some observations which are related to the problem (1). 
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Remark 1. If(x*,y*,z*) €SNSWID(S,%'>B,A,K),by(2)wehavethat 

x*=x*- gi (x*) + [gi (y*) - riATi {ul )] . (3) 
provided K C (//). 

Consequently if 5 is a Lipschitz mapping such that x* e F{S), then it follows from (3) 
that 

X* = Six*) = S (x* -g^ix*) [gi(y*) - nNi {ut,vl)]) . (4) 

By virtue of (4) and Nadler's Theorem (Nadler 1969), we suggest the following iterative 
algorithm. 

Algorithm 1 Let e„ be a sequence of nonnegative real number with e„ ^ 0 as « ^ oo. 
Let r\,r2,rz be three given positive real numbers in (0, 1). For arbitrary chosen initial 
xq e H, compute the sequences {Xn}, {yn} and {z„} in H, such that 

g-iiZn) = Q [g-iiXn) - rsNsiUn,^, V„,3)] , 

gliyn) = [gliZn) - r2N2{Un,2, V„,2)] , Vk > 1 (5) 
Xn+l = (1 - Ol„)X„ + a„S {x„ - gliXn) +Q \gliyn) - ^Afl (««.!. l^w.l)]) ' 

where 

Mn,3 e A^iXn), M„_i,3 & A^iXn-i) : ||m„,3 - M„-l,3|| < (1 + €„)HiAzix„),Aiix„-{)), 

Vh,3 e 53(^h), v„-i,3 e B^{x„-i) : ||v„,3 - v„_i,3|| < (1 + €„)HiBzix„),Bzix„-x)), 

Un,2 e A2iz„),U„-i,2 S A2iz„-i) : \\U„,2 - M„_i,2|| < (1 + €„)HiA2iz„),A2iz„-l)), 

(6) 

Vh,2 e S2(Z„), V„_i,2 e B2iz„-l) : ||V„,2 - V„_i,2|| < (1 + €„)HiB2iz„),B2iz„-l)), 

u„,i eAiiy„),u„-i,i eAiiy„-i) : \\u„,i -m„-i,iII < ii + e„)'H(Aiiy„),Aiiy„-i)), 
Vh,i e Biiy„),v„-i,i e ^iCVn-i) : ||v„,i - v„_i,i|| < (1 + e„)^(Si ()/„), Si 0/„_i)), 

and {«„} is a sequence in (0, 1) and S -.H ^ Hisz mapping. 

Theorem 1. LetK be a nonempty closed and convex subset of a real Hilbert space H and 
<j) -.H ^ (—00, +oo) be a proper convex lower semicontinuous function. LetAi -.H^l^ 
be a m-Lipschitz continuous mapping with m < 1 and Bj : H ^ 2^ be a ai-Lipschitz 
continuous mapping with cr,- < I, i = 1,2,3. Let Ni : H x H ^ H be a pi-Lipschitz 
continuous with respect to the first variable and r]i-Lipschitz continuous with respect to 
the second variable and AT, be Ai-strongly monotone with constant ui > 0 and Bi-relaxed 
monotone with constant > Q,i = 1, 2, 3. Letgi : H ^ H bea Xi-strongly monotone and 
Yi-Lipschitz continuous mapping, i = 1, 2, 3. Let S : H ^ H be a x-Lipschitz continuous 
mapping with Q < t < 1. If SNSWID(S,%^, A,K) D FiS) ^ &, and the following 
conditions are satisfied: 

(i) 



hie 



where 



hi = yi - 2Xi + Yh ^ = 1' 2. 3; 
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(ii) 



|„ ".-ft I ^ ^(vi-Hi)'^-(PilJ-i + m''i)^('ihi)(l-h{) ;_-,90. 



(III) for each i = 1,2,3 

^n,Niiri) + hi ^ ^Njird + hj 
I -hi - I -hi 

where 



< 1, 



(7) 



^„Mn) = ^1 - 2ri(vi - 1;) + rfiipiiii + »7/aO(l + e„))2; 
where M = sup„>i e„. 
(Vvj {a„] c (0, 1) SMC^ that Yl'^o = 

Then the sequences {x„],{y„],{z„],{u„^i],{v„j] suggested by Algorithm 1 converge 
strongly to x*,y*,z*, u*,v* i = 1,2,3 respectively, and (x*,y*,z^,u*,v*) e 
SNSWIDiS, a, ^, A,K), X* e f (5). 

Proof. Let {x*,f,z',i4,vf) e SNSWIDiE,%^,A,K) zndx* e f (5). By (2) and (4) we 
have 

gsiz*) = [g^ix*) - rsNs («*, V*) ] , 

g2(y*) = r; [g2iz*) - r2N2 {ul, V*) ] , (8) 
x* = {l- a„)x*+a„S {x*-gx{x*)+Jl' [gi{y*)-riNx K,v^)]) 

Consequently, by (5) and (6), we have 

\\Xn+l-X*\\ 

= II (1 - an)Xn + a„S (x„ - gi (x„) + [gi (yn) - riNi («„,i, v„,i)]^ - x* \\ 

< (1 -a„)||^„ -;«*!!+ an\\s(^n -gl(Xn)+Q \gliyn) " nNi{Un,l,Vn,l)^ 

- S {f -giix*) +/; [gr{y*) - riNi v^)]) || 

< (1 - oc„)\\x„ -x*\\+ a„x ^\x„ -X* - (giixn) -^i(x*))|| 

+ \\r^'lgl(yn)-nNi(Un,l,Vn,i)]-r^'[gi(y*)-nNi {ut,vt)]\\] 

< (1 II +a«r [||x„ -X* - (giiXn) -^i(^*))|| 

+ Wyn-y* - (gi(yn) -gi(y*m + Wyn-y* - n (iVi(M„,i,v„,i) -m [u^.v*,)) \\] . 

Since Ni(-,-) is pi-Lipschitz continuous with respect to the first variable and 
/71-Lipschitz continuous with respect to the second variable, and Ai is yLti-Lipschitz 
continuous, and Si is cri-Lipschitz continuous, we have 

\\Niiu„,i,v„,i) - Ni Vi) II < pi||m„,i - M^ll + m\\vn,i - vlW 

< pi(l + e„)n(Ai(yn).Ai(y*)) + + e„)mBi(yn).Bi(y*)) 

< Pi/ii(l + e„)||j„-/|| +77iCTi(l + e„)||j„-/|| 

< (Pifii + VK^iXi + e„)\\y„ - y*\\. (10) 
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Since A^i is^i-strongly monotone with constant ui > 0 and Si-relaxed monotone with 
constant > 0, it follows from (10) that 

= \\yn-y*f - 2n (a^i(m«,i. - -^^i {ut.vt) ,yn -/) 

+ rf||Ari(M„,i.v„,i)-A^i(«i.vt) t 
= Wyn-y* f - 2n{Niiu„,uv„,i) - Niiul,v„,i),y„ - y*) 

-2 ri {Ni {u*^, Vn,i) - Ni {ul, v^) , j„ - /) + rl\\Ni(u„,i, v„,i) - Ni {u^, v^) f 

< II J« -/!!'- 2nvi Wyn -y*f+ 2riti \\y„ -y*f 



+ rl{(pilMi + riiaiKl + €„)f\\y„-y*f 
< (1 - 2rii;i+2ri§i + rjUpiHi + ??iori)(l + e„)f) \\y 

n y 



*||2 



i.e., 

hn-y* - n{Ni{Un.i,v„,^)-miu\,v\))f < (ri))2||j„ - / ||l (11) 

where 

^n^in) ■■= y/i- 2ri(ui - li) + r?((piMi + '?i<ri)(l + e„))2. 

Note that 

Il7« - /II = llj« - / - [g2(yn) -g2(y*)] + [gliyn) -g2(y*)]\\ 

< \\yn-y*-[g2(yn) -g2(y*)] II + \\g2(yn) -g2(y*)\\. (12) 

Since ^2 is A.2-strongly monotone and ya-Lipschitz continuous mapping, we have 

II yn-y*-lg2(yn)-g2(y*)]f 

= Wyn -y*W^-2{g2(yn) -g2(y*).yn-y*) + Wg2(yn) -g2(y*)W^ 

< Wyn -y*W^ - 2A2llj„ -/||2 + viwyn -y*W^ 

< (i-2k2 + yi)Wyn-y*W^ 

= (h2fWyn-y*W\ (13) 

where = ^1 - 2A,2 + kI • 
On the other hand, by (2) and (5), we have 

II g2(yn)-g2(y*)W 

= Wl2lg2(Zn) - r2N2(u„,2. V„,2)] [g2(z*) - r2N2 v|)] || 

(14) 

< Wg2{Zn) -g2iz*) - r2iN2iu„,2,V„,2) - N2(ul,vl))W 

< WZn-Z* - (g2(Zn) -^2(z*))ll + WZn - Z* - r2 {N2(Un,2, Vn,2) - !^2 (u^.V^)) ||. 

In view of the assumptions of N2,A2,B2, g2 and by using the same method as given in 
the proofs in (11) and (13), we can obtain that 

||Z„ - Z* - r2(iV2(M„,2, V„,2) - Ar2(M2. V|))||^ < (<^n^M)fWZn - Z*W\ (15) 
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where 

(^nMr2) = ^1 - 2r2(U2 - I2) + r^((P2M2 + '?2fr2)(l + e„))2 

and 

I|Z„ -Z* - te(Z„) -^2(Z*))||2 < ih2f\\Z„-Z*f. 

From (15), (16) and (14), we have 

\\g2(yn) -g2(y*)\\ < i^nMr2)+h2)\\Z„-Z*\\. 

Combining (12), (13) and (17) we obtained 

\\yn-y*\\ < h2\\y„ -/II + {^nN^{r2)+h2)\\z„-Z*\\. 

Observe that 

l|Z«-Z*|| = \\Zn- Z*-[g3{Zn) - g3{Z*)]+[g?.(Zn) - g3{Z*)]\\ 

< \\Zn - Z*-[g^{,Zn)-g^{.Z*)] \\ + ||^3(z») " ^3(z*)l|. 

and in view of (2) and (5), we have 

Wg-iiZn) -g-iiZ*)\\ < \\Xn-X*-[gri{Xn)-g-i{x*)] \\ 

+ \\Xn-X*-r-i {N3(u„,3,V„_3) - N3 (m3,V3))|| 

By using the assumptions on N'iiA^, B3 and ^3, we have 

\\X„ -X* -n (A/3(m„,3, V„,3) - A^3 (M3, V3)) f < (<1)„^3 (rs))^ \\x„ - X* 

where 

^n,N3ir3) = yi - 2r3(U3 - t3) + rl(ip3fJ,3 + r;3(T3)(l + Ch))^ 
ll^«-^* - [g3(Xn) -^3(^*)] 11^ < {h3f\\Xn-X*f. 
\\Zn-Z*-[g3(Z„)-g3{Z*)] f < ih3f\\Z„-Z*f. 

Substituting (21) and (22) into (20), we have 

WgsiZn) -g3iz*)\\ < (<l>«,]V3(r3) + ^3)ll^» -^*ll- 

Combining (19), (23) and (24), it yields that 

I|Z« - Z*ll < ^3||Z« - Z*ll + i^nMn) + h)\\x„ - X*\\. 
This imply that 

WZr.-Z^S ^^-f'l + '^^ lx^-xX 

Substituting (26) into (18) we have 

\\yn -/II < h2\\y„-y*\\ + ''^"'"^^^^^^^^^^I'^-^-'^^^+^^^ lx^ -x1|. 



that is 



... ^ i^nM2(n)+h2X^nMr3)+h3) _ 
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From (11) and (28), we get 

Wyn-y* - ri[Ari(M„,i,v„,i) - A^l(^^^,v^)] II 

(1 -^2)(1 -^3) 

On the other hand, since gi is A. 1 -strongly monotone and yi-Lipschitz continuous 
mapping, we have 

\\X„-X*- (gi(Xn) -glix*))f = \\Xn-X*\\^ + \\gl(Xn) 
-2{Xn-X*,gi(Xn)-gi(x*)) 
< (1 - 2Ai + yI) \\x„-X*\\^ = hl\\Xn-X*\\\ 

i.e., 

Il^« - (gliXn) -^i(^*))ll < hi\\x„-x*\\. (30) 

Similarly, we have 

Wyn-y* - (giiyn) -gi(y*m < hiWyn -y*\\. (3i) 

Substituting (28) into (31), we have 

\\yn-y*-(gi(yn)-gi(y*m 

(1 - /Z2)(l - hs) 

Set 

. ^ (^nMr2) + /i2)(^K,Af3(>'3) + hs) 

Substituting (30), (31), (32) and (33) into (9), we get 

ll*«+i < (1 -a«(l - t(hi+hiln + <^nMri)i„m\Xn-x*\\. (34) 

Since 



< ^1 - 2ri(vi - |„) + rjiipiiMi + ruaiKl +M))2 := ^mM), 

letting I := ^'^'^a ^''(i-fe)(i'-l3)^^"'"^^^ > then we have £„ < I. Therefore from (34) we have 
that 

II -^*ll < (1 - «„(1 - xQi^+h^t + ^NM)f^mxn-x*\\. (35) 
By condition (iii) 



3 



n^^^^<i. (36) 

this imply that 

i < ^ (37) 

that is 

^:=hi + hit + Ojvi (ri)£ < 1. (38) 
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Put 

/7.. = — V* II 

(39) 



tn = 0!„(1 — t3). 

By the assumption that 0 < t < 1, it follows that 

t3 e (0, 1). 

This imply that i„ e (0, 1). From assumption (iv) we have 

00 

^ i„ = 00. 

n=0 

These show that all conditions in Lemma 2 are satisfied. Hence x„ ^ x* n ^ oo. 
Consequently from (26) and (28), we have z„ z* and 3/„ ^ as h ^ oo, respec- 
tively. Moreover since A, is /Lt(-Lipschitz continuous and Bi is cr,-Lipschitz continuous with 
IXi < 1, CT/ < 1, we can also prove that {m„,(} and {v„,(}, j = 1, 2, 3 are Cauchy sequences. 
Thus there exists u*,v* e H such that u„^i — u*,Vn,i — v*, (z = 1,2,3) as « ^ oo. 
Moreover by using the continuity of mappings Ai, Bi,gi, Nitf^, i = 1, 2, 3, it follows from 
(5) that 

g5(z*) = [gzix*) - nNs {ul, vl)] , 

x*=S (x* -gi(x*) +7;^ [gi(y*) - nNi {ut,vl)]j . 

Hence from Lemma 2 it foUows that {x*,y*,z*,u*,v*) e SNSWID(S,a,<B, A,7C). 
Finally we prove that u* eAi(y*) and v* e Bi(y*) Indeed we have 

d(ul,Ai(y*)) = inmul -w\\:w€ Ai(y*)} 

< \\ul - u„,i\\ + d(u„,i,Ai(y*)) 

< \\ut-u„,i\\+mAi(yn).Ai(y*)) 

< - M„,i|| -l-/ii|l:yn -/II ^ Oas« ^ 00. 

That is d(ul,Ai(y*)} = 0. Since Ai(j*) e CB{H), we must have ul e Ai(y*). Similarly 
we can show that e A2{z*),u\ e A-i{x*),v\ e Si()/*),V2 e -82(2*) and Vj e B'iix*). 
This complete the proof. ■ 
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